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t PREFACE 

i 

Work  on  flutter  prediction  has  always  been  based  on  the  modal  representation  of  the 
structure,  which  implies  a linear  model  of  the  aircraft.  A number  of  good  results  have  been 
obtained,  for  many  prototypes,  using  this  method.  Unfortunately,  more  and  more 
difTiculties  have  appeared  during  the  last  few  years,  both  for  aircraft  carrying  large  stores  and 
for  light  aircraft,  where  non-linear  phenomena  make  it  hazardous  to  use  a linear  approach. 

This  pilot  paper,  given  by  Dr  E.Breitbach  to  the  Sub-Committee  on  Aeroelasticity  of  the 
Structures  and  Materials  Panel,  presents  a possibility  for  taking  into  account  some  non- 
linearities  of  the  structure  and  their  effect  on  flutter.  It  was  the  opinion  of  the  Sub- 
committee that  it  should  be  published  by  AGARD  as  an  important  milestone  in  the  way  of 
accurate  flutter  prediction  in  complicated  conditions. 


G.COUPRY 

Chairman.  Sub-Committee 
on  Aeroelasticity 
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EFFECTS  OF  STRUCTURAL  NON-LINEARITIES  ON  AIRCRAFT 
VIBRATION  AND  FLUTTER 


I 


by 

E.  Breitbach 

Deutsche  Forschungs-  und  Versuchsanstalt  fQr  Luft-  und  Raumfahrt  E.  V. 
- Aerodynamlsche  Versuchsanstalt  Gottingen  - 
Institut  fUr  Aeroelastlk 

BunsenstraOe  10,  3400  Gottingen,  Fed.  Rep.  Germany 


SUMMARY 

Experience  has  shown  that  aircraft  structures  are  generally  affected  by  structural  non-linearities. 
The  purpose  of  this  paper  is  to  find  out  the  physical  sources  of  ^e  various  types  of  non-linearities,  to  in- 
vestigate their  Influence  especially  on  the  different  parts  of  the  flutter  clearance  process  and  to  deal  with 
those  methods  which  permit  quantitative  solutions  of  non-linear  aeroelastic  problems. 
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AUi 


time 

structure  point 
frequency 
normal  frequency 

normal  deflection  of  the  quarter -chord  point  of  a wing  section 
rotation  of  a wing  section 
aileron  rotation 

characteristic  angles  of  a non-linear  force  deflection  diagram 

critical  angle  of  limited  amplitude  flutter 

deflection 

normal  mode  deflection 

linear  stiffness  of  an  aileron  hinge 

characteristic  stiffnesses  of  a non-Unear  force  deflection  diagram 

force  or  moment,  F ■ F/C 

generalized  mass 

generalized  stiffness 

gensralized  damping 

generalized  external  force 

damping  loss  angle 

logarithmic  decrement  * 

generalised  coordinate 

linear  stlffbess  energy 


non-linear  stiffness  energy 


nl 


Subscripts 


r.  s 


integers  denoting  the  mode  number 


1.  INTRODUCTION 

As  is  well  known,  aeroelastic  investigations  are  usually  performed  on  the  basis  of  the  simplified  as- 
sumption of  structural  linearity.  However,  it  becomes  more  and  more  evident  that  many  dynamics  and 
flutter  problems  can  only  be  solved  satisfactorily  by  taking  into  account  structural  non-linearities. 

Therefore,  primarily  a survey  is  made  of  the  various  types  of  non-linear  effects  which  may  arise 
on  aircraft  structures.  Beyond  that,  it  is  attempted  to  understand  their  physical  origin  and  to  classify 
them  as  to  their  influence  on  ground  and  flight  test  methods  as  well  as  on  the  analytical  flutter  predic- 
tion. Emphasis  is  placed  on  methods  which  are  suitable  for  treating  non-linear  vibration  and  flutter 
problems  with  special  regard  to  systems  with  several  degrees  of  freedom  formulated  in  terms  of  calcu- 
lated or  measured  modal  characteristics. 

The  applicability  of  such  methods  is  Illustrated  by  a number  of  experimental  and  analytical  results. 


2.  SURVEY  OP  RELEVANT  STRUCTURAL  NON-LINEARITIES  ON  AIRCRAFT 

Non-linearities  as  arising  in  aircraft  structures  can  generally  be  subdivided  into  distributed  non- 
linearities  which  are  continuously  activated  throughout  the  whole  structure  by  elasto -dynamic  deforma- 
tions and  into  concentrated  ones  which  act  locally  lumped  especially  in  control  mechanisms  or  in  the 
connecting  parts  between  wing  and  external  stores.  Some  characteristic  types  of  these  two  classes  of 
non-linearities  are  described  and  discussed  in  the  following  sections. 


2.  1 CONCENTRATED  NON-LINEARITIES 


The  control  mechanisms  of  hand  operated  aircraft  are  affected  by  different  types  of  strong  concen- 
trated non-linearities  resulting  from 

• back-lash  in  the  linkage  elements  of  the  control  system 

• solid  friction  in  control  cable  and  push  rod  ducts  as  well  as  in  the  hinge  bearings 

• kinematic  limitation  of  the  cor‘rol  surface  deflection 

• application  of  spring  tab  systems  provided  for  relieving  pilot  operation. 

These  effects  generally  occurring  in  combined  forms  can  beet  be  demonstrated  by  means  of  results 
measured  on  real  aircraft  structures.  Fig.  1 shows  I 

the  aileron  hinge  moment  of  a glider  versus  hinge 
angle  in  case  of  static  moments  symmetrically  act- 
ing in  the  sense  opposite  to  the  regular  operation  of 
the  aileron  system.  The  result  represents  the  sim- 
plest form  of  a low  damped  back-lash.  The  anti- 
symmetrical  force  deflection  curve  resulting  from 
static  loading  in  the  regular  operation  sense  of  the 
aileron  system  is  given  in  Fig.  2.  It  is  a hysteresis 
type  diagram  with  an  elastic  and  nearly  frictionless 
slope  for  forces  below  the  critical  slip-stick  point 
where  static  friction  changes  to  sliding  friction.  The 
latter  causes  an  hysteresis  Increasing  up  to  the  ma- 
ximum aileron  deflection,  beyond  which  the  force 
deflection  curve  is  characterized  by  nearly  friction- 
less  elastic  deformation.  Resulting  from  this  non- 
linearity the  resonance  frequency  of  the  antlsymms- 
trlcal  aileron  vibration  decreases  from  about  7 Hz 
at  low  amplitudes  up  to  a minimum  of  2.5  Hz  at 
about  15  degrees,  s.  Fig,  3,  The  resonance  frequen- 
cy Increases  again  along  the  dashed  line  for  ampli- 
tudes exceeding  the  maximum  aileron  stroke.  The 
force  deflection  diagram  of  a glider  rudder  in  Fig.  4 
is  almost  exclusively  characterized  by  solid  friction 
whereas  elastic  deformations  are  negligible. 

Quite  similar  non-linear  effects  due  to  back- 
lash and  solid  friction  can  be  observed  on  aircraft 
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Fig.  1;  Aileron  hinge  moment  of  a glider 
versus  hinge  angle,  symmetrical 
loading 


Fig.  2:  Aileron  hinge  moment  of  a gli- 
der versus  hinge  angle,  antl- 
symmetrlcal  loading 


Fig.  3:  Resonance  frequency  of  the  anti- 
symmetrical  aileron  vibration 
versus  hinge  angle 


with  external  wing  stores.  In  [1]  and  [2]  special  emphasis  Is 
placed  on  the  Investigation  of  the  single  point  suspended  un- 
der wing  stores  of  a combat  aircraft  with  sweepable  wing. 
Main  purpose  of  these  studies  Is  to  find  out  the  limits  of  line- 
arisation and  their  possible  consequences  to  test  and  analy- 
sis. 

Another  type  of  a concentrated  non-llnearlty  Is  occurring 
In  spring  tab  systems  deliberately  Introduced  Into  the  control 
mechanism  In  order  to  relieve  pilot  operation.  Both  the  me- 
chanism of  a spring  tab  system  and  the  appertaining  force 
deflection  curve  are  sketched  In  Fig.  5 and  Fig.  6 respective- 
ly, assuming  the  hinge  stiffness  Fq  of  the  main  control  sur- 
face to  be  zero.  The  stiffness  of  the  spring  tab  system  Is  re- 
latively high  In  a rather  small  range  around  the  origin  of  the 
diagram  resulting  from  a special  pre -tension  of  the  tab 
spring  Ff . At  a certain  tab  deflection  the  pre-tension  cea- 
ses and  the  stiffness  suddenly  drops  to  a much  lower  value. 
The  force  deflection  diagram  of  the  complete  control  me- 
chanism Is  a combination  of  the  tab  spring  F-p  and  the  stiff- 
ness Fc  between  control  stick  and  control  surface  which  Is 
generally  characterized  by  back-lash  and  solid  friction.  As 
outlined  In  one  of  the  following  sections,  such  spring  tab  sys- 
tems have  already  proved  to  be  very  susceptible  to  a danger- 
ous kind  of  divergent  flutter. 

Power  operated  controls  as  schematically  shown  In  Fig.  7 
are  also  affected  with  non-linear  effects.  This  might  be  il- 
lustrated for  Instance  by  the  resonance  frequency  of  the 
F 104  G aileron  system  plotted  in  Fig.  8 as  a function  of  the 
aileron  deflection,  [3].  Generally,  It  cam  be  assumed  that 
the  dynamic  behaviour  of  power  control  mechanisms  Is  cha- 
racterized by  a combination  of  different  sorts  of  non-linear 
effects  such  as  outlined  In  the  following; 

s As  pointed  out  In  [3]  amd  [4]  the  complex  highly  damped 
stiffness  of  ttie  hydraulic  actuator  Itself  Is  a non-linear 
function  of  the  vibration  amplitude  as  well  as  of  the  actua- 
tor preload  amd  of  the  position  of  the  Jack  piston. 

• The  linkage  mechanism  between  the  actuator  body  amd  the 
atlrcraft  structure  Is  aiffected  by  considerable  non-linear 
effects  due  to  solid  friction  and  back-lash  especially  In  the 
hinge  bearings. 

• In  case  that  vibratory  motions  of  the  atlrcraft  structure  re- 
sult In  relative  motions  between  piston  and  cylinder  of  the 
servo  valve,  actuator  forces  are  induced  which  also  vary 
with  the  vibration  atmpUtude.  This  kind  of  dynamic  inter- 
action can  lead  to  considerable  stability  problems. 


Fig.  4;  Rudder  hinge  moment  of  a glider 
vereus  hinge  angle 


In  addition  to  all  that,  the  dynamic  stiffness  of  hydraulic 
amtuators  Is  usually  a frequency  dependent  function  to  be  ta- 
ken Into  account  as  a very  Important  part  of  the  flutter  clea- 
rance process,  see  [5]. 


2. 2 DISTRIBUTED  NON-UNEARITIES 

The  vibration  behaviour  of  aircraft  structures  Is  influ- 
enced not  only  by  the  concentrated  non-linear  effects  des- 
cribed above  but  also  by  the  so-called  distributed  ones  which 
are  Induced  by  elastic  deformations  in  riveted,  screwed  and 
bolted  connections  as  well  as  within  the  structural  compo- 
nents themeelvee.  Because  of  the  great  number  of  rivets, 
screws  and  bolts  It  can  be  assumed  that  the  resulting  damp- 
ing and  stiffness  non-Unearltiee  are  more  or  less  continuous- 
ly distributed  throughout  the  structure.  In  consequence  of  this 
special  property,  the  effect  of  these  distributed  non-Unearl- 
tles  on  the  normal  mode  shapes  can  moetly  be  regarded  as 
negligible. 

Experience  has  shown  that  the  normal  frequencies  are 
weakly  decreasing  functions  of  the  vibration  amplitude.  Quite 
contrary  the  overall  damping  values  can  undergo  much  hitter 
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variations  against  amplitude.  Two  typical  results 
of  a ground  vibration  test,  carried  out  on  the  air- 
craft F 104  G,  [6],  are  plotted  in  Fig.  9, 

As  found  out  in  special  damping  measurements 
on  a fiber  reinforced  composite  box  spar,  [7], 
structural  damping  of  modern  design  components 
depends  not  only  on  the  vibration  amplitude  but  also 
on  frequency.  Fig.  10  shows  the  results  of  the  fun- 
damental bending  normal  mode,  the  normal  fre- 
quency of  which  was  modified  from  235  up  to  88  Hz 
by  attaching  additional  masses.  Thus,  the  stress 
distribution  over  the  test  structure  could  be  kept 
invariant. 


3.  INFLUENCE  OF  STRUCTURAL,  NON-LINEARI- 
TIES ON  TEST  AND  ANALYSIS 

The  survey  above  demonstrates  that  there  are  a 
lot  of  different  types  of  structural  non-linearities  to 
be  taken  into  account  in  ground  vibration  test  as 
well  as  in  flutter  analysis  and  flight  flutter  test. 

At  first  sight,  the  above  described  distributed 
non-linearities  do  not  seem  to  be  of  great  impor- 
tance for  aeroelastic  investigations.  However,  many 
difficulties  as  occurring  particularly  during  ground 
vibration  and  flight  flutter  tests  can  be  attributed  to 
these  effects. 

Thus,  it  has  been  shown  in  (8)  that  weak  distri- 
buted stiffness  non-linearities  in  concurrence  with 
equally  weak  manufacture -conditioned  structural  a- 
symmetries  can  be  identified  as  physical  sources  of 
the  well-known  phenomenon  of  amplitude  dependent 
normal  mode  asymmetries. 

Another  very  interesting  problem  concerning 
the  measurement  of  generalized  masses  is  investi- 
gated in  (9|.  The  stucfy  comes  to  the  conclusion  that 
measured  generalized  masses  can  be  considerably 
affected  by  measuring  errors  due  to  small  distri- 
buted stiffness  non-linearities.  Suitable  means  to 
overcome  this  problem  are  proposed. 

As  to  the  Influence  of  distributed  damping  non- 
linearities,  it  is  worth  mentioning  that  certain 
changements  of  the  critical  flutter  speed  resulting 
from  non-linear  damping  effects  can  be  expected 
only  in  cases  of  mild  flutter.  However,  damping  be- 
comes far  more  Important  in  case  of  many  dynamic 
response  problems  because  the  dynamic  amplifica- 
tion of  an  externally  excited  structure  is  inversely 
proportional  to  the  structural  damping. 

Regarding  the  so-called  concentrated  non-line- 
arities it  is  self-evident  that  in  most  cases  simpli- 
fied linearized  approaches  must  lead  to  inadmis- 
sible errors  and  misinterpretations.  To  settle  these 
problems  some  propositions  have  been  made  espe- 
cially in  view  of  a more  sophisticated  preparation 
and  execution  of  the  ground  vibration  test. 

According  to  [10]  non-linear  dry  friction  effects 
on  control  surfaces  of  hand  controlled  aircraft  can 
be  largely  reduced  by  the  application  of  an  auxiliary 
periodic  high  frequency  excitation.  In  this  way  an 
aircraft  structure  can  be  artificially  linearized  for 
frequencies  much  lower  than  that  of  the  auxiliary 
excitation. 

If  the  modal  synthesis  approach  is  utilised  the 
detrimental  influence  of  concentrated  non-lineari- 
ties on  the  accuracy  of  ground  vibration  tests  can  be 


Fig.  5:  Sketch  of  a spring  tab  mechanism 


Fig.  6;  Schematic  force  deflection  dia- 
gram of  a pre -tensed  spring  tab 


Fig.  7;  Sketch  of  a power  operated  con- 
trol surface 
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Fig.  8;  Resonance  frequency  of  the  F 104  G 
aileron  vibration  versus  trailing 
edge  amplitude 
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eliminated  by  means  of  three  different  procedures,  which  guarantee  the  measurement  of  a set  of  ortho- 
gonal normal  modes: 


• During  the  ground  vibration  test  the  mechanism  of  the  control  surfaces  can  be  linearized  by  re- 
placing the  non-linear  elements  by  linear  artificial  ones.  As  shown  in  [11],  the  real  non-linear 
effects  can  be  subsequently  reintroduced  into  the  flutter  analysis  by  means  of  modal  superposi- 
tion. 

• In  accordance  with  [12]  during  the  test  the  control  surfaces  can  be  fixedly  clamped  to  the  adja- 
cent wing  or  tail  structure  in  such  a way  that  relative  deflections  are  suppressed.  This  implies 
the  elimination  of  the  control  surface  degrees  of  freedom.  The  real  vibration  behaviour  can  be 
completely  described  in  the  subsequent  analysis  by  introducing  the  control  surface  degrees  of 
freedom  in  addition  to  the  measured  ones.  The  non-linear  effects  of  the  control  mechanisms  can 
be  taken  into  account  using  an  analytical  approach  quite  similar  to  the  one  described  in  [11], 

• As  described  in  [13]  the  aircraft  structure  can  be  investigated  in  the  ground  vibration  test  with- 
out control  surfaces.  In  this  way  non-linear  effects  cannot  be  activated.  The  normal  modes  of 
the  control  surfaces  are  separately  measured.  The  mathematical  model  of  the  complete  struc- 
ture can  be  subsequently  established  by  means  of  a modal  coupling  approach  taking  into  account 
the  non-linear  properties  of  the  coupling  mechanism. 


Fig.  9:  Structural  damping  and  resonance  fre- 
quency of  a typical  F104G  normal 
mode  versus  reference  amplitude 


Fig.  10;  Structural  damping  of  a fiber  reinforced 
composite  box  spar  aa  function  of  vibra- 
tion amplitude  and  frequency 


Every  one  of  the  three  procedures  quoted  above 
requires  the  force  deflection  diagrams  of  the  con- 
trol surface  mechanisms.  It  has  proved  to  be  the 
most  effective  and  accurate  way  to  determine  these 
properties  experimentally. 

As  to  flight  flutter  test,  concentrated  non-line- 
arities must  be  taken  into  consideration  as  inherent 
properties  of  the  test  object.  Mamipulations  as  ap- 
plied to  ground  vibration  measurements  in  order  to 
facilitate  them  are  not  possible  here. 

The  methods  to  interprets  and  evaluate  flight 
flutter  test  results  reliably  can  be  best  improved  by 
means  of  detailed  analytical  and  wind  tunnel  investi- 
gations on  systems  with  realistic  non-linearities. 

In  the  following  the  governing  aeroelastic  equa- 
tions of  an  aircraft  affected  by  concentrated  non-li- 
nearities will  be  formulated  in  terms  of  modal  cha- 
racteristics. In  order  to  give  an  impression  of  the 
influence  of  some  typical  non-linearities  on  the  flut- 
ter stability  the  results  of  an  analogue  compaiter 
flutter  simulation  will  be  presented.  Furthermore, 
the  applicability  of  a slightly  modified  numerical 
approach  based  on  the  principle  of  "harmonic  ba- 
lance" will  be  checked  up  by  means  of  the  results  of 
a wind  tunnel  flutter  test  on  a non-linear  wing  ailer- 
on model. 


3. 1 FORMULATION  OF  THE  NON-LINEAR  AERO- 
ELASTIC EQUATIONS 


Provided  that  the  ground  vibration  test  is  car- 
ried out  on  an  artificially  linearized  aircraft  struc- 
ture the  governing  equations  can  be  written  on  the 
basis  of  the  measured  modal  data  as  follows 


with 


^r^  °rB^l  * W • «r 

r.s  ■ 1.2.. ...N 


(1) 


generalized  masses 
generalized  stiffnesses 
generalized  damping  coefficients 

q^  generalized  coordinates,  4r(^)  • <le>>ote 

the  first  and  second  derivation  with  respect 
to  time 
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Q generalized  unsteady  aerodynamic  forces  depending  on  ^ . Mach  number  and 

*"  Reynolds  number. 


The  generalized  stiffness  can  be  expressed  by 

K » 4ir^f^M  (2) 

r r r 

with  the  normal  frequency  fr  of  the  r-th  normal  mode.  The  solution  q^  of  Eq.  {!)  can  be  retransformed 
into  the  corresponding  geometrical  displacements  u(P,t) , which  are  defined  as  a function  of  time  t and 
of  the  structure  points  P,  by  means  of  the  series  expansion  of  the  measured  normal  mode  functions 
Ur(P) 

u(P.t)  • 2 u JP)  q.(t)  : s-1,2 N . (3) 

(8)  ® 

As  mentioned  before,  in  the  ground  vibraUon  test  the  non-linear  stiffness  elements  characterized  by  the 
force  deflecUon  functions  FmOm)  are  replaced  by  artificial  linear  stiffnesses  Cn . Subscript  p indi- 
cates the  control  surfaces  taken  into  consideration  and  denotes  the  hinge  angle  of  the  (»-th  control 
surface. 

The  equations  of  motion  of  the  non-linear  system  can  be  established  by  adding  to  Eq.  (1)  the  La- 
grange term 

^ (AU  , - AU.)  : r*l,2....,N  (4) 

dq  nl  I 


with  the  stiffness  energies  AUni  and  AUi  stored  in  the  non-linear  and  the  linear  springs.  This  opera- 
tion is  equivalent  to  the  reestablishment  of  the  real  non-linear  conditions  by  analytical  means.  The  en- 
ergy terms  AUj  and  AUni  are  defined  as  follows: 


AU  • i S C 

^ Mm)  " 


AU  » I / FO  )dd  : M = 1.2 L 

(M)  0 


(5) 


Applying  to  a series  expansion  similar  to  that  formulated  in  Eq.  (3)  results  in 


^ ^us‘*s  • ® “ 1.2,. ...N  (6) 

M (s) 

where  stands  for  the  hinge  angle  of  the  M-lh  control  surface  as  measured  in  the  s-th  normal  mode 

of  the  linearized  system. 

Insertion  of  Eq.  (6)  into  Eq.  (5)  leads  to 


(7) 


and  with 


to 


^ B dq 
(.)  “ 


s • 1,2, 


,N 


(8) 


fi., 


AU 


nl 


(M)  0 


s • 1,2, 


..,L 

..,N 


(9) 


If  the  Lagrange  operation  is  applied  in  accordance  with  Eq.  (4)  to  Eq.  (7)  and  (9)  and  the  result  is  added 
to  Eq.  (1)  this  leads  to  the  aeroelastlc  equations  of  motion  of  the  non-linear  aircraft 


M 4 4^ 
r r 


(■) 


D 4 

re  a 


•f  K q ♦ 
r^r 


r.s  • 1,2,....N 
ft  ■ 1, 2, ....  L . 


(10) 


Eq,  (10)  is  based  on  the  plausible  hypothesU  that  the  dynamic  displacements  of  the  non-linear  system  can 
be  formuUted  in  the  same  way  as  those  of  the  linear  system,  namely  as  a superposition  of  the  normal 
modes  of  the  linearized  test  configuration.  This  formulation  is  very  advantageous  because  not  only  the 
modal  daU  , Drg  and  K,  btrt  also  the  unsteady  aerodynamic  forces  Qr  con  be  taken  over  from  the 
linear  system  unchanged. 

The  equations  can  be  solved 

e In  the  domain  by  mesne  of  an  analogue  computer  or  by  numerical  step  by  step  integration. 

e In  the  frequency  domain  by  means  of  iteration  procedures. 


The  derivation  of  the  governing  Eqs.  (10),  especially  with  regard  to  the  unsteady  aerodynamic  forces, 
is  dealt  with  more  in  detail  In  (11). 


7 


3,2  FLUTTER  SIMULATION  OF  A SIMPLE  NON-LINEAR  WING  AILERON  SECTION 

In  order  to  investigate  the  Influence  of  control  surface  non-linearities  on  the  flutter  stability.  In  [11] 
a real  time  flutter  simulation  by  means  of  an  analogue  computer  was  carried  out  for  a simple  wing 

aileron  section  with  a span  width  of  100  cm.  The 
unsteady  aerodynamic  forces  were  determined  by 
means  of  the  incompressible  strip  theory  including 
Wagner's  function  for  non-uniform  motions.  As 
sketched  in  Fig.  11  the  system  is  characterized  by 
the  following  degrees  of  freedom 

wing  normal  deflection  h 
wing  torsion  o 

aileron  rotation  . 

The  modal  characteristics  Mr,  fr  and  the  normal 
mode  shapes  hr,  Or , fir  "linearized"  basic 

configuration  are  summarized  in  Table  1. 

The  aeroelastic  stability  curves  resulting  from 
three  different  types  of  non-linear  force  deflection 
characteristics  are  plottet  in  Fig.  12  to  14.  The 
stability  boundaries  are  presented  in  form  of  the 
ratio  of  the  critical  aileron  hinge  angle  to  a reference  angle  versus  flight  speed.  By  way  of 
comparison,  the  flutter  boundary  of  the  linear  basic  system  is  demonstrated  in  form  of  a vertical  line  at 
V~  300  km/h. 


Normal 

Mode 

fr 

(Hz) 

Mr 

2, 

[kg  cm  ] 

hr 

[cm] 

"r 

[rad] 

Pt 

[rad] 

r a 1 

3.58 

10.  18 

1 

0.  66  10'® 

0.50  10*^ 

2 

15.  33 

26.01 

1 

-0.039 

-0.21 

3 

18.  13 

116.72 

1 

-0.045 

1.  17 

Fig.  11:  Sketch  of  a wing  aileron  section  with 
three  degrees  of  freedom 


Table  1:  Modal  characteristics  of  the  "linearized"  basic  configuration  of 
the  wing  aileron  section 


As  shown  in  Pig.  12  non-linear  flutter  resulting  from  a back-lash  type  hinge  stiffness  leads  to  a li- 
mited amplitude  flutter  which  is  not  explosive  but  may  cause  a long-term  process  of  material  fatigue. 


VElXITYlkm/hJ 

Fig.  12;  Effect  of  back-lash  on  the  flutter 
stability  of  a wing  aileron  section 


Incase  that  the  force  deflection  curve  is  cha- 
racterized by  a combination  of  back-lash  and  solid 
friction  the  flutter  boundary  line  is  shifted  to  con- 
siderably higher  speed  values,  see  Fig.  13.  This 
stabilizing  effect  is  characteristic  for  all  hystere- 
sis type  hinge  stiffnesses. 

In  contrast  to  that,  flutter  resulting  from  de- 
creasiiig  or  softening  non-linear  springs  can  be 
much  more  dangerous.  It  can  be  seen  from  Fig.  1 4, 
that  this  spring  tab  type  hinge  stiffness  is  very  in- 
sidious, because  flutter  does  not  start  until  a cer- 
tain critical  amplitude  limit  is  exceeded,  for  in- 
stance in  consequence  of  a gust  excitation.  Above 
this  critical  limit  the  vibrations  are  steadily  in- 
creasing until  destruction  of  the  aircraft. 

As  described  in  [14]  in  a qualitative  manner  for 
a great  number  of  examples,  the  flutter  behaviour  of 
systems  with  non-linear  control  surface  mecha- 
nisms is  highly  dependent  on  the  manner  the  criti- 
cal flutter  speed  of  the  appertaining  "linearized* 
systems  varies  against  changes  of  the  linear  hinge 
stiffness  as  shown  in  Fig.  15  for  the  above  investi- 
gated wing  aileron  section.  Functions  of  that  kind 
can  help  towards  a first  qualitative  asessment  of 
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the  flutter  stability  of  non-linear  systems. 


3.  3 FLUTTER  ANALYSIS  AND  WIND  TUNNEL  TESTS  ON  A NON-LINEAR  WING  AILERON  MODEL 

Recently,  in  the  DFVLR/AVA  low  speed  wind  tunnel  flutter  tests  were  carried  out  on  a half  span 
wing  aileron  model  with  a non-linear  aileron  hinge  stiffness.  The  main  dimensions  of  the  model  can  be 
seen  from  Fig.  16.  In  accordance  with  the  above  described  procedure  the  ground  vibration  test  was  per- 
formed on  a basic  configuration  with  a "linearized"  aileron  hinge  stiffness.  The  measured  modal  cha- 
racteristics M]-,  fr.  the  damping  loss  angles  yr  and  the  appertaining  normal  mode  shapes  are  sum- 
marized in  Fig.  17.  The  relation  between  the  logarithmic  decrement  Or  > decaying  oscillation  and 
the  loss  angle  yr  ie  • 


Fig.  13;  Effect  of  hysteresis  on  the  flutter 
stability  of  a wing  aileron  section 
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Fig.  14:  Effect  of  a softening  spring  tab  type 
stiffness  on  the  flutter  stability  of  a 
wing  aileron  section 


Fig.  15;  Flutter  speed  of  the  "linearized*  wing  aileron 
section  as  function  of  changes  of  the  linear 
hinge  stlffhess 


§ 


% 


50cm 


W7cm 


Fig.  16;  Main  dimenaions  of 

wind  ttinnel  test  model 


fi -5.08  Hz 

Hj-802kgcin^ 

7j-0.0i7 
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f,  -.3.33Hz 
M,  -0.63  kg  cm^ 
7,  =0.135 


fj  = 7.42Hz 
Mj=3.2kgcm^ 
73  = 0.025 


Fig.  17:  Normal  mode  characteristics  of  a wing  aileron 
wind  tunnel  test  model 


Fig.  IB;  NoC'linear  force  deflection  diagram  of  the  aileron  of  the 
wind  tunnel  test  model 


The  force  deflection  curve  of  the  wind  tunnel  test  configuration  shown  in  Fig.  18  is  characterized  by 
backslash  of  t 3 mm  in  combination  with  a small  hysteresis. 

Fig.  19  shows  the  flutter  boundary  of  the  non>linear  wind  tunnel  teat  model  in  comparison  with  cal- 
culated results  in  form  of  the  critical  amplitude  ratio  aileron  as  a four  valued  function  of 

velocity  V.  The  two  full  lines  indicate  so-called  limited  amplitude  flutter,  whilst  the  dashed  lines  des- 
cribe indifferent  states  only  to  be  determined  by  calculation.  Thus,  along  the  upper  dashed  flutter  bound- 
ary even  extremely  small  external  disturbances  result  in  increasing  or  decreasing  vibrations  converging 
against  one  of  the  full  lines. 

According  to  the  strength  of  external  disturbances  due  to  gust  or  manoevre  loads  the  limited  ampli- 
tude flutter  turns  up  along  the  upper  or  the  lower  full  line. 

During  the  tests  it  turned  out  that  the  damping  values  were  highly  dependent  on  the  momentary  vi- 
bration amplitude  so  that  a reliable  damping  measurement  was  Impossible.  This  non-linear  effect  can  be 
observed  in  flutter  tests  on  wind  tunnel  models  as  well  as  in  flight  flutter  testa,  irraspectlve  of  the  damp- 


w 
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Ing  measurement  method  applied. 

Regarding  the  risk  of  life  and  safety  of  an  aeroelastic  system  resulting  from  limited  amplitude  flut- 
ter Fig.  19  demonstrates  that  the  flutter  amplitude  ratio  varies  from  ^ic/^a~  V • 12  m/s  to 

^k/^a~  ® finally  reaches  the  asymptotic  case  **  V ■ 19  m/s  predominant- 

ly defined  by  the  stiffness  C2  beyond  the  flat  hysteresis  part  of  the  force  deflection  curve  in  Fig.  18. 


FLUJJER  SPEED  [m/$] 


Fig.  19:  Comparison  of  the  measured  and  the  calculated  Fig.  20:  Sketch  of  a hysteresis  type 

flutter  boundary  of  the  non-linear  wind  tunnel  force  deflection  diagram 

test  model 


The  calculation  of  the  non-linear  flutter  boundary  was  performed  on  the  basis  of  the  principle  of  the 
"harmonic  balance"  by  means  of  a rather  simple  numerical  approach.  In  accordance  with  [IS]  the 
elastomechanic  properties  of  a non-linear  force  deflection  diagram  can  be  approximately  described  by 
equivalent  linear  coefficients,  an  equivalent  stiffness  coefficient  CeO)  end  an  equivalent  damping  loss 
engle  yeiP) . which  can  be  calculated  from 

2» 

CeO)  ■ ^ y F(/Jcosv),  - ^ u sin«>)  cos^  d^ 
to  ‘0 

(12) 

2w 

" w'fc  (it)  f u •Inip)  alaifi  dtp  . 

® (P-0 

The  integration  variable  of  these  amplitude  dependent  functions  is  defined  by  * ut  with  w ■ 2xf . 

The  solution  of  the  integrals  for  a hysteresis  type  force  deflection  diagram  like  that  sketched  in 
Fig.  20  yields 

[2FQ(sln<PQ- cos^pj)  - * costp^)  ♦ 

* f *®2'  sin2v>j+  2<Pj,.  2<pj)  ♦ f 

’ (IS) 

- 1(02-01)  (cos2v)g*cos29i)] 


tPjj  ■ arc  cos  -j- 


tpj  • arc  sin  -j- 


The  equivalent  aileron  hinge  stiffness  of  the  non-linear  hinge  mechanism  of  the  wind  tunnel  model  is 
plotted  in  Fig.  2 1 as  a function  of  the  amplitude  ratio  /3/0a*  ^ further  important  result  indispensable 
to  the  determination  of  the  flutter  boundary  of  the  non-linear  system  is  given  in  Fig.  22  illustrating  the 
variation  of  the  flutter  speed  against  hinge  stiffness,  which  is  assumed  to  be  linear  for  this  calculation. 

Finally,  the  flutter  boundary  of  the  non-linear  wing  aileron  system  as  shown  in  Fig.  19  could  be  de- 
termined from  a comparative  evaluation  of  Fig.  21  and  Fig.  22. 

Calculated  and  measured  results  agree  very  well  and  thus  demonstrate  that  the  application  of  the 
principle  of  the  "harmonic  balance”  is  very  efficient  even  in  case  of  hl^ly  non-linear  systems. 

The  "harmonic  balauice*  approach  can  also  be  applied  to  systems  with  more  than  one  concentrated 
non-linearity.  Starting  from  the  general  equations  of  motion  (10)  this  problem  can  be  solved  by  numeri- 
cal iteration. 


AMPLIIUDERAriO  ft/p. 

Fig.  21:  Variation  of  the  equivalent  aileron  hinge 
stiffness  of  the  wind  tunnel  model  ver- 
sus hinge  angle  ratio 


Fig.  22;  Flutter  speed  of  the  "linearised*  wing 
aileron  wind  tunnel  model  as  function  of 
changes  of  the  linear  hinge  stiffness 


4.  CONCLUSION 

Both  a survey  and  a classification  of  structural 
non-linearities  as  occurring  on  aircraft  have  shown 
that  a variety  of  non-linear  effects  have  to  be  taken 
into  account  in  the  course  of  aeroelastic  investiga- 
tions. The  physical  properties  of  distributed  and 
concentrated  non-linearities  have  been  exemplified 
by  means  of  a number  of  typical  results  measured 
on  real  aircraft  structures.  Possibilities  to  over- 
come detrimental  non-linear  effects  on  ground  and 
flight  flutter  tests  have  been  discussed. 

The  influence  of  some  different  types  of  aileron 
hinge  non-linearities  on  the  flutter  stability  has 
been  illustrated  by  the  results  of  a flutter  simula- 
tion carried  out  on  a simple  wing  aileron  section  in 
the  time  domain.  Furthermore  the  flutter  behaviour 
of  a wing  aileron  wind  tunnel  model  with  a non-li- 
near aileron  hinge  mechanism  has  been  investigated. 
The  test  results  have  been  successfully  compared 
with  calculated  ones  obtained  by  a numerical  proce- 
dure on  the  basis  of  the  principle  of  the  "harmonic 
balance". 

As  to  future  investigations  emphasis  should  be 
placed  on  the  following  problems: 

a Improvement  of  the  ground  vibration  test  concept 
for  instance  by  artificial  linearisation  of  the  test 
structure  and  by  subsequent  introduction  of  the 
non-linearities  into  the  aeroelastic  equations. 

a Investigation  of  the  non-Unear  and  frequency  de- 
pendent (tynamic  properties  of  power  operated 
controls:  incorporation  into  the  aeroelastic  equa- 
tions. 

a Investigation  of  weak  distributed  non-linearities 
as  to  their  effect  on  ground  vibration  test  (with 
special  regard  to  modern  experimental-numeri- 
cal methods)  and  flight  flutter  test. 

a Amendment  of  flight  flutter  test  methods  aiming 
at  a more  reliable  detection  of  the  non-linear 
flutter  boundaries. 

a Further  development  of  numerical  methods  to 
solve  non-linear  flutter  problems  with  several 
intensively  interacting  concentrated  non-lineari- 
ties. 

a Elaboration  of  suitable  methods  to  predict  unstea- 
dy aerodynamic  forces  due  to  arbitrary  motions 
for  real  wings  with  finite  aspect  ratios  and  for  all 
flight  speed  ranges. 
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14.  Abstract 


mork  on  flutter  prediction  has  always  been  based  on  the  modal  representation  of  the  structure, 
which  implies  a linear  model  of  the  aircraft.  A number  of  good  results  have  been  obtained,  for 
manv  prototypes,  using  this  method.  Unfortunately,  more  and  more  difficulties  have  appeared 
during  the  last  few  years,  both  for  aircraft  carrying  large  stores  and  for  light  aircraft,  where 
non-linear  phenomena  made  it  hazardous  to  use  a linear  approach.  This  paper  presents  a 
possibility  for  taking  into  account  some  non-linearities  of  the  structure  and  their  effect  on  flutter. 
It  is  considered  to  be  an  important  milestone  in  accurate  flutter  prediction  in  complicated 
conditions,  yjy 

Paper  rpesented  at  the  45th  Structures  and  Materials  Panel  Meeting.  Voss,  Norway,  September 
1977. 


